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Abstract 

We consider the dynamics of charge carriers in single-layer graphene that are subject to random 
temporal fluctuations of their mass gap. The optical conductivity is calculated by incorporating 
the quantum-stochastic time evolution into the standard linear-response (Kubo) theory. We find 
that, for an intermediate range of frequencies below the average gap size, electron transport is 
enhanced by fluctuations. At the same time, in the limit of high as well as low frequencies, 
the conductivity is suppressed as the variance of gap fluctuations increases. In particular, the dc 
conductivity is always suppressed by a random temporal mass with nonvanishing mean value and 
vanishes in the zero-temperature limit. Our results are complementary to those obtained recently 
for static random-gap disorder in finite-size systems. 

Keywords: Graphene, Temporal mass-gap fluctuations. Quantum-stochastic time evolution. 
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1. Introduction 

Graphene is a newly accessible nanomaterial lUl S Bl^ consisting of a single sheet of 
carbon atoms forming a two-dimensional honeycomb lattice. The valence (highest occupied) 
and conduction (lowest unoccupied) bands of graphene are touching at the K, K' = -K (and 
equivalent) points in the Brillouin zone and exhibit a conical shape of their dispersion in the 
vicinity of these special (Dirac) points H 0, [sd . Graphene is a promising candidate for ap- 
plications in future micro- and nanoelectronics due to its excellent mechanical characteristics, 
scalability to nanometer sizes, and the ability to sustain huge electric currents [9, 10]. Recent ex- 
periments have confirmed that the charge carriers in graphene indeed behave like massless Dirac 
fermions f^Jv^. 

Although graphene exhibits excellent conducting properties, the absence of an energy gap 
poses a challenge for realizing conventional semiconductor device operations in this material. 
A possible way to alleviate this situation is to induce a gap in the electronic spectrum at the 
Fermi energy by breaking the discrete sublattice symmetry of the honeycomb structure. In recent 
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experiments lll2Lll3ll was reported that hydrogenation of the graphene sheet breaks the sublattice 
degeneracy. Reversibility of this process enables switching between conducting and insulating 
regimes of a graphene sample, and spatial addressability should make the fabrication of hybrid 
conducting-insulating graphene samples possible. 

The opening of a uniform gap destroys graphene's metallic state, as the ac (optical) conduc- 
tivity vanishes for photon energies smaller than the band gap. In real samples, the gap may be 
random, e.g., because of structural variability in the physical system. For example, the bonds with 
hydrogen atoms could undergo temporal and spatial variations. To understand this practical issue, 
and also for more fundamental reasons, studies of static random-gap disorder in Dirac-fermion 
systems have been performed ifTil [l5 , 16, 17], where a staggered potential was considered as 



a possible model for breaking the sublattice symmetry of the honeycomb lattice. Our work 
presented here augments these previous investigations by considering the effect of a spatially 
uniform mass gap that fluctuates in time. Besides clarifying the respective nature of temporal 
and spatial randomness in the mass gap, which is an issue of basic interest, the scenario consid- 
ered by us would be reahzed in graphene samples where sublattice-potential-inducing agents are 
undergoing large-scale random variations. Below we argue that a more general spatio-temporal 
random variation of a staggered sublattice potential in graphene can be represented, to leading 
order, by a Markovian model where the low-energy Dirac-fermion theory has a random mass 
gap described by white-noise fluctuations. We use this model to derive our stochastic equation 
of motion for the disordered-graphene system and incorporate the nontrivial time evolution into 
the framework of linear-response theory. Thus a generalized Kubo formula is obtained for the 
optical conductivity, enabling us to unambiguously determine how electric transport is aff'ected 
by the mean value and the variance of the fluctuating potential. This is a useful result because 
the optical conductivity is an experimentally accessible quantity |21L |22|] that depends 

strongly on the electronic properties of the material. In particular, it can be used to determine the 
size of a band gap 123,1 since transport is suppressed for photon energies smaller than the gap. 
We find that the dc conductivity at finite temperature is monotonously decreasing as a function 
of the variance of mass-gap fluctuations. However, for an intermediate range of frequency below 
the cut-off' equivalent to the average gap, fluctuations enhance the conductivity. This observa- 
tion should be contrasted with the related increase of the conductivity with variance of static 



mass-gap fluctuations in graphene samples of finite size II15LI17I1 



The remainder of this paper is organized as follows. Our theoretical method is outlined in the 
following Sec. |2] where a generalized linear-response formula for the conductivity of a system 
subject to random temporal fluctuations is derived. Results obtained from application of this 
formalism to randomly gapped graphene are presented in Sec. [3] We explore dependencies on 
various experimentally controllable parameters and discuss the relation of our results to previous 
work. Our conclusions are given in Sec. 



2. Theoretical method 

2.1. Fluctuating gap model 

Quantum (and classical) systems experience dissipation and fluctuations through interaction 
with a reservoir or environment llist . Here, the hydrogenization of the graphene sample induces 
states, which we consider as a reservoir The interaction with the states of the reservoir gener- 
ates fluctuation, and a short explanation of the technical details is given in [Appendix A[ In the 
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context of the open quantum system (subsystem and reservoir), we consider the tight-binding 
Hamihonian for electrons in graphene subject to a sublattice-staggered potential, 

•H = -f 2 + H.C.) + 2 Vi,a a'lh^ + J] Vtb bl,ks, (1) 

(ij).s i,s i,s 

where aj^, a,,, annihilates (creates) an electron with spin s (s =T> i) on site / on sublattice A 
(an equivalent definition is used for , on sublattice B), t ^ 2.8eV is the nearest-neighbor 
hopping energy between different sublattices. V,,,, and Vi^t are spin-independent potentials with 
Vi^a - mi on sublattice A and V/j, - -m, on sublattice B. These potentials break the sublattice 
symmetry of the single-layer graphene. 

As a consequence of the interaction with the reservoir, m, is a time-dependent random vari- 
able at lattice site /, having a mean value ((m,(r))) = '« and a variance {{(mi(t)-m){mj(t')-m)}} - 
g^5ij6(t - f). The double expectation value is just the average over the local probability distribu- 
tion function related to the lattice site /, followed by averaging over the whole lattice probability 
distribution function. In a more general case [24], values for m,(f) at different lattice sites could 
be correlated. We make the assumption that the random variables are uncorrected for different 
lattice sites and the value of the average gap (averaging locally) is the same for each lattice site. 
This random-mass variable can be modeled as m,(f) = m + g^i(t), where ^,(f) represents temporal 
white-noise fluctuation at lattice site /, and {{^i{t)^j{t'))) - Sijd(t - t'). 

We remind the reader that the static disorder models discuss only spatial correlations be- 
tween different lattice sites. The aforementioned description given by us could deal with local 
(temporal fluctuations) and spatial correlations, too. However, in the present paper we consider 
a complementary case to the static disorder models. Our model for randomness is given by the 
relation m,(f) = m + g^i{t), and it is clear that an average over the local fluctuations leads to a 
spatially uncorrected variable m. 

The usual derivation 101 of the low-energy continuum model for quasiparticles in graphene 
yields the Hamiltonian 

H(k) = hv(cr,k, + o-,,k,.) + m{t)cr.ik, (2) 

with multiplication operators k^, kv, Ik (that become numbers k^, ky, 1 in the plane-wave repre- 
sentation), and cr, denoting the Pauli matrices for the pseudo-spin degree of freedom, v is the 
Fermi velocity, which has a value ^ 10^ m/s. The m{t) represents the randomness introduced in 
the tight-binding model; it is equal to m + g^it), where ^(f) is the usual white noise. The white 
noise, ^(f), dependence on the ^j{t) is defined by Yij^ji^)^'^^'^^^^' - ^(0<^(k' - k), where k, k' 
are the variables of the Fourier transformation, and Rj is the position vector of lattice site j. 
We assume disorder potentials to be weak (inducing small momentum scattering) such that the 
inequivalent Dirac points associated with the K and K' valleys in graphene's band structure are 
uncoupled f23\ and, thus, contribute additively to electronic transport. This allows us to consider 
the presence of randomness for each valley separately. The fact that electrons in graphene also 
carry a real spin introduces an additional double degeneracy of all eigenvalues for the Hamil- 
tonian (|2|. In the following, we study the dynamics due to the single-valley, spin-less model 
defined by Eq. (|2]) and absorb degeneracy factors into our basic unit of conductivity. 

2.2. Quantum-stochastic time evolution for the density matrix 

As the Hamiltonian (|2|i contains a random variable, the time evolution induced by it cannot be 
expressed by the usual equations. Instead, a quantum-stochastic calculus has to applied. A brief 
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introduction to that is presented in the following. l|26 , 27 1 An isolated quantum system evolves 
in a unitary fashion. A physical quantity that is given at time f = by an observable A, will be 
described at time f > by A(t) = U^{t)AU{t), where U(t) is a unitary operator for each point in 
time t. The unitary operator is generated by the Schrodinger equation: 

dUit) i . 

= --Hom(t), (3) 

where the (time dependent) Hamiltonian H^)(t) is a self-adjoint operator for each t. The Schrodinger 
equation tells us that in a short time interval dt, the unitary operator changes at time t like 

dU(t) = -^Ho(t)U(t)dt. (4) 

n 

The Hamiltonian in Eq. (|2]) can be written in the form of H{k) - Ho(k) + Hi^(t), where 

H^ik) = hv(o-_^k^ + CTykj.) + mo-,lh, (5a) 
Hi = gcrX (5b) 

The white noise ^(f) is the formal derivative of a Wiener process Wt (continuous everywhere but 
difFerentiable nowhere), a Gaussian random variable with zero mean value 

M(W,) = 0, (6) 

and variance t 

MiWf) - M(W,f = t. (7) 

The stochastic calculus applied here will be based on the Wiener process Wt. The above proper- 
ties are represented in the differential equations as 

M(dW,) = 0, d^W, = dt, d"W, = 0, n > 2, (8) 

all of which will be applied in further calculations. Consider a stochastic processes X, governed 
by a stochastic differential equations 

dX, = V,dt + D^dW,, 

where physical meaning of Vx is the drift and of is the diffusion. The Ito rule for a smooth 
function f{x) states: 

df{X,) = f'{X,)dX,+ ^-f"{X,){dXtf, 

df(x) „ d^f(x) 
/ W = f"{x) = -f^, X 6 R, 

dx dx-^ 

(dX,f = Dldt. 

The quantum Ito rule fl^ is the same as the Ito rule, only the noncommutativity of the operators 
is considered in addition. 

Now, we apply the properties of the stochastic calculus, and we see that the unitary operator 
changes within an infinitesimal short time interval dt at time t as follows: 
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dU(t) = ( --Ho(k)dt - —Hldt - -HxdW, I U{t). (9) 



The equation of motion of the density matrix is given by p(f) - U(t)pU'(t), and a straight- 
forward calculation yields 

dpit) = -^[Ho(k),pit)]dt-^[Hu [Hup(t)]]dt-^[Hup(t)]dW,. (10) 

dr 

The current operator in our system is defined by 7^ — e-^, where the time evolution of the 
coordinate operator satisfies 



dr^it) = -[i/o(k),r/f)]fl?f- — [^^i,[^^i,?;.(0]]fl?f+ T['^i'^A'(0]fl?W'f (H) 



is a derivative operator in the k-space and, in the sublattice representation, is given by i2®;g|-. 
This operator commutes with Hi, which has a structure of cr^®!];. These conditions and Eq. ( fTTT i 
define the current operator as 7'^ - ^[Ho(k), r^] - f Due to the operator structure in 

Eq. (|2]i, the single-particle eigenstates \n) of this model can be written as a direct product of a 
plane wave in configuration space with a spinor: |n) = |k) (g) |+)k (^,^2^. Here + labels the 
electron and hole bands, respectively, and the spinor wave function depends on wave vector k. 
The Cartesian current components are y\ = evcr^ and y, = every ISOll . The definition of the 
equilibrium density matrix in the spinor space is pol±)k = /(fk,±)|c)k, and i/ol±)k = fk,±l±)k, 
where / is the Fermi-Dirac distribution function and ek,± = ± V(^^)^|kP + frfi- This equilibrium 
state is the average of all equilibrium realizations defined by po - e^l^^^i-) iYj-(e^^^^^^). 

2.3. Generalization of linear-response formalism 

We employ the linear-response (Kubo) formalism llsill and divide the system's Hamiltonian 
into the part H{k), which governs the evolution in Eq. (fTOl i. and 6H, the perturbation associated 
with an external electric field E. For simplicity, we take the latter to be constant in space and 
assume the field to be applied between t = -00 and f = 0. The perturbation Hamiltonian is 
6H - -eE ■ f e"^'. The equation of motion of the system with the added external field is: 

dp(t) = -kHo(k),p(t)]dt - ksH,p(t)]dt -:^[Hu [Hupmdt - kHup(t)]dW, 

£pit)dt - UsH,pit)]dt - UHup(t)]dW,. 

h h 

Within linear-response theory, we can linearize p - pQ + 6p, where po is the system's equilibrium 
density matrix. Keeping only linear terms in Eq. (fTSl l. we get 

ddp = JLSpdt - U5H,pQ\dt - [Hi,pQ]\dt - UHuPo + 5p]dW„ (12) 

with using [i/o(k),po] = as well as [6H, 6p] - 0. Introducing Ap = e^-^'6p yields 

dAp = e-^'{-^[5H,po] - ^[Hi, [Hi,po]]jflff - ^^-^'{[HuPo + e^'Ap])dW,. (13) 

Note that Ap and 6p have the same value at f = 0, and both vanish at f = -00. Integration yields 

Spit = 0) = J dt e-^'l-^[SH,po] - ^[Hu [i/i,po]]j+lim J dW, (e"^' (^[i/i,po]) - ^[Hu 

(14) 



The exponential factor in Eq. (fl4l l ensures convergence of the time integral, making it unneces- 
sary to introduce the phenomenological adiabatic damping parameter employed in conventional 
linear-response theory llsil . 



2.4. Generalized conductivity formula 

We use (fT4l i to calculate the double expectation value (over the basis of the one-particle 
Hilbert space and over the ensemble of realizations for the stochastic process) for the current 
operator j and use the property M(t/W,) - 0. We remind the reader that our stochastic variable 
models the local gap fluctuations. An average over this randomness yields the mean conduc- 
tivity 1 32] that is typically calculated and is our quantity of interest here. Dividing the current 
expectation value by |E| yields the optical conductivity 

%J — OO 

with the kernels 

= -^Tr{]^e-^'([er„po])} , (16) 

L^v(f) = Tr||ie-^'J-^[^,,[^i,po]])| • (17) 

The Tr symbol stands for taking the trace over 2x2 matrices and performing the integration 
J c/^k. As the integrand in Eq. (fTTl i is an odd function of both and A:,., the integration over the 
polar angle yields zero and, hence, the quantity L^v(f) vanishes. 

We remind the reader that the current flowing through the system was not affected by the fluc- 
tuations, see Eq. (fTTT l and the discussion afterwards. Otherwise, the variance of the fluctuations 
must be added to the conductivity formula. This constant would shift the universal conductance 
of the graphene sheet, and would refute the experimental results. This is not the situation here, 
which shows the correctness of our model. 

The calculation of K^y is straightforward, and using the Laplace transform to solve for the 
dynamics, we find 



h J I 



where Res stands for the sum of residues of the integrand. The conductivities related to the 
kernels and Kyy are identical, and we present the effect of the fluctuating gap for these 
particular conductivities. We have to mention that the Hall conductivities cr „ and cr,, ^ are equal to 
zero, even in the presence of the fluctuations. Without an external magnetic field, the contribution 
of the two valleys of the honeycomb-lattice band- structure cancel each other 

Introducing polar coordinates for k and performing the angular integration yields 

e^ P |k| d\k\ ^ f N(Mm,g,z) ,,\ 

K^y(t) = —\ — Res^ _ -e''\. (19) 

ft Jo 47r [D{\k\,m,g,z) J 

For simplicity, the Fermi velocity v has been absorbed into |k|, and we find 

dm , \ 



N{\k\,m,g,z) = |k|2 Vot2 + |k|2 (W + 4\k\^ + z^ - ATz) x (" 



de '^=^^- 



+ z(z-4r)(2m2 + |krt[/(ek,-)-/(ek,+)] , (20) 
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Figure 1 : Optical conductivity of single-layer graphene with a weakly fluctuating mass gap. Results shown are obtained 
for fixed chemical potential fi = IQ hgT and gap-fluctuation variance SF = 0.1 kgT. The different curves correspond to 
rh/ksT = 5 (solid), m/lcBT = 10 (dashed), and mlkgT = 15 (dotted). The value of is marked by a horizontal line. 
As can be seen, the ac conductivity exhibits a jump at twice the average band-gap value in when the latter is greater than 
or equal to the chemical potential. 



D{\\i\,m,g,z) = (m^ + |k|2)5 [\6zT^ - 8(2|k|2 + z^)T + z(\vr? + \\\^ + z^)] , (21) 

where the parameters F = g^/2?i^ and m - fhlh have been introduced. The cubic factor in the 
denominator dTTT i has three roots z,, which give the poles in eq. ( fT9] l. 

3. Optical conductivity of gapped graphene 

For small F, the roots are to lowest order zi = and Z2,3 = ±2; + nf-. Using this and 
performing the time-integration in the limit F ^ 0, the /nfra-band contribution 



Ji I lis''' 

= -5(oj) \ 
and the inter-hand contribution 



= 77 T, z — ~ 2m), (23) 

o-o 8 aP- cosh(j^)-Hcosh(5|^) 



to the conductivity of the uniformly gapped single-layer graphene are found 112311 . The Dirac- 
delta peak in the /nfra-band conductivity is due to elastic transitions, which are only possible at 
finite temperature and/or when the chemical potential is bigger than the average gap. The scale 
factor (To - 4e^/h accounts for spin and valley degeneracy. 

In the following, the conductivity is calculated numerically for finite values of F and T from 
Eq. ( fTSl l with Eq. (fT9] l. We use ksT as our unit of energy. Figs. [T] and |2] show the ac (optical) 
conductivity, whereas Figs.[3]|4]and|5]show results for the dc conductivity. 

The interplay between average-gap size and chemical potential exhibits two regimes, as can 
be seen by studying the uniform gapped model [result given by Eq. (l23T l1 and also from Fig. [1] 
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Figure 2: Optical conductivity of single-layer graphene with a fluctuating gap. Top panel: the chemical potential is fixed 
at the neutrality point (jilksT = 0), and the average gap is m/irgr = 5. Bottom panel: fi = WkgT and m = SksT. The 
different curves coiTespond to hT/ksT = 0.1 (solid), KT/kBT = 1 (dashed), and hF/ksT = 5 (dotted). The value of 
7t/8 is marked by the horizontal line. As the value of F is decreased, the curves resemble more closely the result found 
for uniformly gapped graphene. Note the intermediate range of low frequencies where an increase in randomness (i.e., 
variance of fluctuations) results in an enhanced conductivity. 



The uniform gap model has a vanishing ac conductivity for photon energies ha> < 2fh, and a 
jump at hcij = 2m. The height of the jump depends on the chemical potential and the average 
gap. If the chemical potential is much bigger than the average gap then we have a small jump 
height whereas, in the opposite case, the jump is more noticeable. The ability to determine the 
average gap size by measuring where the jump of the optical conductivity occurs depends also 
on the variance of the fluctuations. In any case, setting the chemical potential as low as possible 
would be favorable for the detection of the gap-induced jump in the ac conductivity. 

Figure |2] shows the ac conductivity for a number of different variances of gap fluctuations for 
two values of the chemical potential fi: tuned to the neutrality point (top panel) and for a large 
value of jj (bottom panel). While the situation with a small variance of gap fluctuations is quite 
similar to the result found for a system with a uniform gap, the increase of fluctuations generates 
sizable conduction in the frequency range w < 2m. For both values of the chemical potential, 
a region emerges where an increase in fluctuations results in an increased conductivity. Thus 
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Figure 3: The dc conductivity of randomly gapped grapliene plotted as a function of gap variance T. The average gap 
is fixed at m/kgT = 2.5. The different curves are for ii/ksT = (solid) and fj/ksT = 5 (dashed). A decrease in the dc 
conductivity is found for increasing F. This decrease is slower when the chemical potential is bigger than the average 
gap- 



it appears that larger fluctuations will facilitate electric transport in randomly gapped graphene 



sheets. A similar result has been reported for static random-gap disorder in graphene fl5l 11711. 
However, in our case, this behavior occurs only within a limited range of finite frequencies. The 
parametric dependence on the fluctuation strength is reversed at low frequencies, in particular 
also for the dc conductivity. For high frequencies, the conductivity approached the universal 
value 7r/8. The detailed shape of the saturation depends on both F and m, with higher values 
pushing convergence to higher frequencies. 

The dc conductivity's dependence on the parameters F, m, and the chemical potential is 
shown by Figs.[3]|4]and|5] Increasing F leads to a decrease in the dc conductivity, which means 
that the F simulate the same effect as life-time broadening due to inelastic scattering. The rate 
of the decrease depends on the relative magnitudes of fi and m, as can be seen in Fig. [3] If the 
chemical potential is bigger than the average gap, the dc conductivity tends to zero slower with 
any increase of F. Increasing the average gap also leads to a decrease in the dc conductivity. The 
minimal value of the dc conductivity occurs at the charge neutrality point, fi - 0, and strongly 
depends on the value of F and m. 

Up until now, we considered electron transport at a finite temperature. To discuss the be- 
havior in the limit when T ^ 0, we have to perform this limit already in Eq. (|T9) and redo the 
calculations following after that. In the case of T —> Q and = we are able to compare our 
results to the works by Ziegler OlSll . These works assume that the static disorder average destroys 
the intra-band conductivity, which we also found for the above assumptions. We are able to go 
further and in the case of T ^ andu < m the intra-band conductivity is still found to be zero. 
In the analytical works by Ziegler OlSll . the inter-band conductivity contains a Heaviside function 
and has a contribution at w = 0. The role of the Heaviside function is to separate the insulating 
regime from the metallic one. This is a point where we have a significant difference, because our 
Heaviside function is smeared by the F parameter In the limit w = with in i^Q and fi < m, our 
inter-band conductivity is always zero. 
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Figure 4: The dc conductivity of randomly gapped graphene plotted as a function of average gap for different values of 
gap variance T. The different curves correspond to hT/ksT = 0.5 (solid), HF/kBT = 1 (dashed), and hr/ksT = 1.5 
(dotted). The chemical potential is fixed at the neutrality point. As the average gap size is increased, the dc conductivity 
gets suppressed. 



The role played by the quantity F in the present work is formally similar to a decoherence 
parameter (also denoted by T in a previous work [Issl] by us) measuring the effect of an ever- 
present environment. The parameter F is the coefficient of a double commutator that enters to 
the system's time evolution. In the present case, the double commutator contains the d-~ operator, 
in contrast to our previous work li33t where it involved the o".v and o",, matrices. The physical 
meaning of T is also different: here it is related to the randomness of a mass-gap-inducing degree 
of freedom such as hydrogen atoms, while it was used to describe the properties of a current 
detector previously. We would also like to point out that, in general, the gap fluctuations could be 
temperature-dependent, as most mechanisms for inducing a gap (e.g., by structural modifications 
such as hydrogenation) will be affected by a variation of temperature. As we have introduced 
the fluctuating gap in a phenomenological fashion, our theory would apply to the ensemble of 
temporal fluctuations at a fixed temperature. 

4. Conclusions 

We have calculated the stochastic evolution of the density matrix for charge-carrier dynam- 
ics in single-layer graphene subject to a fluctuating sublattice-staggered potential. We derived 
a generalized Kubo formula to study the effect of temporal mass-gap fluctuations on the con- 
ductivity. The variance of the fluctuations introduces a source of damping and thus makes the 
converged adiabaticity parameter frequently used in Kubo formula calculations superfluous. See 
also related work [33]. Mixing of the intra-hand and inter-hand contributions to the ac and 
dc conductivities strongly affect its parametric dependence on the variance of the fluctuations 

- 2h^r. A system with a uniform gap (i.e., vanishing variance F) is insulating, as seen from 
Eq. ( l23T l. For the more realistic case of a fluctuating gap, a finite range of frequencies is found for 
which an increase in randomness (i.e., variance of gap fluctuations) results in an increase of the 
ac conductivity. See Fig.|2] Such fluctuation-enhanced transport is not observed in the dc limit. 
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Figure 5: The dc conductivity of randomly gapped grapliene. plotted as a function of chemical potential for different 
values of gap variance F. The different curves were obtained for hT IUbT = 0.5 (solid), hTIhgT = 1 (dashed), and 
hr/ksT = 1.5 (dotted). The average gap is fixed at m/ksT = 2.5. The minimal value of the dc conductivity occurs at 
fj = Q and decreases with increasing T. 



where increases in either the average size or the variance of the mass gap lead to a suppression 
of conductivity at any finite temperature. In the zero-temperature limit, our model of randomly 
gapped graphene exhibits insulating behavior when the average gap is finite 1 14, 16] (but, as in 
the previous work Il4il . a finite dc conductivity is found when the average gap vanishes). 

The discussed diff'erences between different random gap models can be found out in the prop- 
erties of the conductivity, the only quantity retained after a measurement. 

The author has profited from helpful discussion with K. Ziegler, A. Sinner and U. Ziilicke. 
This work was supported by a postdoctoral fellowship grant from the Massey University Re- 
search Fund. Additional funding through BMBF project QK_QuOReP is gratefully acknowl- 
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Appendix A. Stochastic dynamics 

A quantum master equation describes the dynamics of an open quantum system (subsys- 
tem-Hreservoir). The Markovian quantum master equation is given by a first-order differential 
equation for the density matrix of the subsystem, 

^Ps(t)^-Cps(t), (A.l) 
at 



where the most general form for the generator X is given by the Lindblad equation 1134 13511 . 

-CPs ^-^[Hs,Ps]+Yj ^* i^^kPsAl - ^AlAkPs - ^psAlAi^ . (A.2) 

The first term of the generator represents the unitary part of the dynamics generated by the 
Hamiltonian .^5 of the subsystem. The operators A^ are the Lindblad operators and the quantities 
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yic have the dimension of an inverse time. The operators and the quantities are derived from 
the dynamics of the total system (subsystem and reservoir) in various approximation schemes. 
The states of the reservoir are traced out during the procedure. 

It is clear that the eq. (lA.ll i cannot be used for the derivation elaborated in the work of 
Semenoff [TJ . The master equation for the subsystem can be reformulated in terms of a stochastic 
process for the subsystem's wave function. This idea is the so-called unravelling of the master 
equation ll36l 371. 

Now, we consider our microscopic model of the total system, 

H = H« +H, + H\ 

(i,j).s i,s i,s 

= Yj'^'^l^'-^^"h.s,a-Yj'^lKshs®Hl...,t^ (A.3) 



H. 



where aj_,, a, j annihilates (creates) an electron with spin s (s -1, X) on site / on sublattice A (an 
equivalent definition is used for b''.^, bj^s on sublattice B). W is the Hamiltonian of the states, 
considered as reservoir Hi is the interaction Hamiltonian and W, and W, . , are functions of 

^ l.i,s.a 1,1, s,b 

creation and annihilation operators related to the states of the reservoir coupled to sublattices A 
and B. The interaction Hamiltonian is chosen such that the sublattice symmetry is broken. A, is 
a spin independent coupling constant, m,- is the strength of the sublattice symmetry breaking on 

lattice site /. 

The usual Bom-Markov master equation [27] can be derived by considering that the time 
scale of the hopping is much more faster than the relaxation of the reservoir states and the cou- 
plings are weak: 

^ = £p(t) = -^[H^ps(f)] -^j^ Tr, [H,it), [Hjis),ps(t)®pr]] ds, (A.4) 

where Trr stands for the partial trace over the reservoir states, pr is the density matrix of the 
reservoir Inserting Eq. (IA.3) into the master equation (Eq. ( IA.4I )) we obtain 

i,s 

H'j = fl;>,,-%^/,„ (A.5) 

where gi is related to the reservoir correlation functions, and we considered that these correlations 
are the same for sublattice A and B, too. 

The unravelling of the above master equation gives the following stochastic Hamiltonian 

H^-tYj {a'lbj,., + H.C.) + ^ mi(t) ci\ a^ - ^ '«/(0 ^- (A-6) 

{i,j),s i,s i,s 

where 

m,(f) = m + gm- (A.7) 



^(f) stands for the standard white noise. 
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